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ABSTRACT 

Cadenced optical imaging surveys in the next decade will be capable of detecting 
time- varying galaxy-scale strong gravitational lenses in large numbers, increasing the 
size of the statistically well-defined samples of multiply-imaged quasars by two orders 
of magnitude, and discovering the first strongly-lensed supernovae. We carry out a 
detailed calculation of the likely yields of several planned surveys, using realistic dis- 
tributions for the lens and source properties and taking magnification bias and image 
configuration detectability into account. We find that upcoming wide-field synoptic 
surveys should detect several thousand lensed quasars. In particular, the Large Syn- 
optic Survey Telescope (LSST) should find more than some 8000 lensed quasars, some 
3000 of which will have well-measured time delays. The LSST should also find some 
130 lensed supernovae during the 10-year survey duration, which is compared with 
~ 15 lensed supernovae predicted to be found by a deep, space-based supernova sur- 
vey done by the Joint Dark Energy Mission (JDEM). We compute the quad fraction in 
each survey, predicting it to be ~ 15% for the lensed quasars and ~ 30% for the lensed 
supernovae. Generating a mock catalogue of around 1500 well-observed double-image 
lenses, as could be derived from the LSST survey, we compute the available precision 
on the Hubble constant and the dark energy equation parameters for the time delay 
distance experiment (assuming priors from Planck): the predicted marginalised 68% 
confidence intervals are <j(wq) = 0.15, a(w a ) — 0.41, and a(h) = 0.017. While this 
is encouraging in the sense that these uncertainties are only 50% larger than those 
predicted for a space-based type-la supernova sample, we show how the dark energy 
figure of merit degrades with decreasing knowledge of the the lens mass distribution. 
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1 INTRODUCTION 

The discovery of st rong gravitational lensing in Q0957+561 
(Walsh ct al. 197§) opened up the possibility of using strong 
lens systems to study cosmology and astrophysics. Roughly 
a few hundred strong lens systems produced by massive 
galaxies have bee n discovered to date, with quasars (e.g 



Inada et al.l |200ST) . radio sources (e.g.. Myers et al 



Browne et aljl2003), and galaxies (e.g . JCabanac et alj 



2003 



2007 



Faure et alj |200S| ; iBolton et ail 120081 ) appearing as source 



populations. These samples of strong lenses have been ex- 
tensively used to constrain dark energy, the density profile 
of lensing gala xies, and t he evolution of massive ellipticals 
(see, e.g.. iKocha nck 2006, for a review). 

It is expected that large new samples of strong lenses 
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will be obtained in future wide-field im aging surveys (e.g., 
iKuhlen et al. l l2004lMarshall. et al.|[2005l l. For instance, the 
Large Synoptic Survey Telescope (LSST; [ivezic et al.ll2008h 
will observe a 20,000 deg 2 region with a final limiting mag- 
nitude of r ~ 27.5, which is considerably wider and deeper 
than existing optical i maging surveys su ch as the Sloan Dig- 
ital Sky Survey CSDSS: lYork et al.|[200oh . More importantly, 
future imaging surveys such as LSST's will pay particular 
attention to the acquisition of time-domain data, mostly in 
order to study transient objects including supernovae (SNe). 

In this paper, we present detailed predictions of the 
numbers of time-variable sources we can expect to be 
strongly lensed in ongoing and future time domain optical 
imaging surveys. We focus on time- variable sources for two 
reasons. First, time domain data enable us to i dentify stron g 
lenses by taking a dvant age of time variability (|Pindoill2005l h 
iKochanek et~afl (|2006l l proposed to find strong lenses by 
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looking for "extended variable sources" in time-domain data, 
and argued that this technique should be very efficient due 
to the low levels of contamination in the difference images. 
Second, time-variable sources allow us to measure time de- 
lays between multiple images; these time delays contain rich 
information on the lens potentials as well as cosmology (e.g 



Witt ct al 
SchechteT 



2000l; iKochanekl 12002]; ISaha fe William; 



- .20051; lO guri 200 7al: ICongdon et al.l I200S 
Keeton fc Moustakasll2009l : lKeetonll2010l ). 



2003 



2010 



The obvious and traditional example of such a strongly 
lensed time- variable source is a quasar (QSO). Of ~ 100 
gravitationally lensed quasars currently known, time delays 
have been measured for only < 20 (see IOgu"rHl20 07a. for a 
recent compilation). A major advance offered by the upcom- 
ing imaging surveys is their potential not only to discover 
new lenses but also monitor them as well. Some of the sur- 
vey telescopes we will consider do have larger static surveys 
planned, but we do not consider them on the grounds that 
if they are wide enough to be competitive, they will be lim- 
ited by the required monitoring follow-up. Synoptic surveys 
(such as those designed to discover supernovae) should pro- 
vide some estimate of the lens time delays - we are careful 
to point out where this may not be the case. 

The scope of several planned surveys allows us also to 
consider SNe as time-variable sources. While a distant SN 
highly magnified by a foreground massive c luster has re- 
centl y been identified |Stanishev et al.l 120091 ; iGoobar et al.l 
120091 ), no multiply -imaged SN has yet been discovered. How- 
ever, the possibility of discovering such strongly-lensed SNe 
by future time-domain surveys such as the Joint Dark En- 
ergy Mission (JDEM) an d the LSST h as been pointed out 
(iPorciani fc Madaul l200d; iHolj [200ll; IGoobar etafl |2002| ; 



lOguri et al.ll2003l ; iMortsell et al.1 120051 ). A particularly en- 
ticing feature of lensed SNe is that, if they are type-la, the 
"standardizable candle" nature of the peak brightness allows 
a direct determination of the lensing magnification factor, 
which breaks the deg eneracy between the len s potential and 
the Hubble constant (|Oguri fc Kawano| [2003). 

Indeed, we expect any large sample of time delay lenses 
to be useful in constraining cosmological parameters. The 
idea of using strong gravitational lens systems with time- 
vari able sources to measure the Hubble constant is an old 
one (jllcfsdal 1964|); in fact, the time delays between images 
provide a way of measuring the "time-delay distance" to the 
lens. This quantity is a combination of the angular diame- 
ter distances to the lens, source, and between the two, and 
while primarily sensitive to Hq, it does also depend on the 
other cosmological parameters (see e.g.. lSuvu et al.|[2010l . for 
a recent high-precision measurement). Attempts to measure 
the Hubble constant from t he statistica l analysis of the cur- 
rent lensed quasar sample l|Saha et al.ll200a lOguril l2007al : 
IColesl 2008) have yielded values cons istent with constraints 
from other cosmo l ogica l probes (e.g.. iFreedman et al.ll2001 



P_ 

Bonamente et al.l 120061 ; Komatsu et al.l 1200 9; Rics s et al.l 



2009). In future, larger samples covering a wide range of 



lens and source redshift permutations should provide inter- 
esting constraints on the dark energy equation of state as 
well as Hubble constant (iLewis fc Ibatal | 2002l ; iLinderl [20041 ; 
iDobke et af1l2009l ; ICoe fc Moustakasll2009l ). In this work, we 
revisit the forecast constraints on dark energy parameters 
for the LSST, using new, realistic predicted distributions of 
strong lenses. 



This paper is organised as follows. In Section [2] we de- 
scribe the ingredients of our lens abundance calculation, in- 
cluding our assumptions about both the lens and source pop- 
ulations. We then present our predictions for various survey 
depths and areas in Section [3] We then present Monte Carlo 
realisations of some fiducial survey lens catalogues in Sec- 
tion [4] and use these to explore the potential of some of 
these in constraining cosmological parameters via the time 
delay distances to the lenses in Section [S] After a brief dis- 
cussion (Section [6]) we conclude in Section [7] We take as 
the parameters of a fiducial cosmological model fi m = 0.26, 
Ode = 0.74, h = 0.72, w = -1 and w a = 0. We denote the 
angular diameter distance from observer to lens as D\, from 
observer to source as D s , and from lens to source as Di s . 



2 PREDICTING THE NUMBERS OF LENSED 
QUASARS AND SUPERNOVAE 

We compute the expected number of strong lensing in fu- 
ture surveys in the usual way, integrating the lensing cross 
section of each galaxy over th e galaxy population and red- 
shift (e.g.. iTurner et al.|[l984l) . As source objects, we con- 
sider both QSOs and SNe, for which measurements of time 
delays between multiple images will be possible due to their 
time-variable nature. 



2.1 The population of lensing galaxies 

An appropriate modeling of the lens galaxy population con- 
stitutes an essential part of the predi ction for len s ing ra tes. 
We mostly follow the model used in I Oguri" et all l|2008l ) to 
compute the expected number of strong lensing by mas- 
sive galaxies in future surveys. Throughout the paper we 
conservatively consider only elliptical galaxies (E/S0) as 
lensing objects; these l ikely comprise ~ 80% of the to- 
tal lensing probability l|Turner et al.l Il984l; iFukugita et al 



19921 ; lKochanek1ll996l ; Ictiac 2003; Oguri 2006; Mollcr et al 



2007). We assume that the mass distribution of each ellip- 
tical galaxy is described by a Singular Isothermal Ellipsoid 
(SIE), whose convergence is given by 



K(x,y) = 



A(e) 



0Ein 



47T 



fv\ 2 As 
\c) Ds' 



(1) 



(2) 



where v is the one-dimensional velocity dispersion of the 
lensing galaxy. It has been shown that this model describe 
the mass distributio ns of observed strong lensing ellipti- 
cal galaxies very well ( Rusin & Ma 2001; Treu fc Koopmansl 
20041; iRusin fc Kocha nck 20051; iKoopmans et al.ll2006l . l2009l ; 
Gavazzi et al.ll2007l ). The parameter A(e), the so-called dy- 
namical normalisation, is related t o the three- dimensional 
shape of lensing galaxies. Following [Cha^ (|2003l ). we assume 
that there is an equal number of oblate and prolate galaxies, 
and adopt the average value of the no rmalisations i n eac h 
of the oblate and prolate cases. As in lOguri et all |2008), 
we assume a Gaussian distribution for the ellipticity e, with 
a mean of 0.3 and dispersion of 0.16. The distribution is 
truncated at e = and e = 0.9. 

In addition to the SIE model lens galaxy, we include 
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external s hear to account f or the effect of the lens environ- 
ment fe.g jKochanek 1991; Witt fe Maoll 19971 ; iKeeton et al.l 
1 19971 ). The lens potential of external shear is given by 



y ) cos 28 1 + -yxy sin 28^ 



(3) 



We assume the magnitude of the external shear, 7, follows 
a log-normal distribution with mean 0.05 and dispersion 
0.2 dex, which is broadly consistent with the level of exter- 
nal shear expected from ray-tracing in TV-body s imulations 
l|Holder fc Schechted 12003. ; foala! fc Watsonl ,20041 ). The ori- 



entation of the external shear, # 7 , is assu med to be random . 
Although the external convergence re.g jQguri et al. 1,20051 ) 
is also vital particularly for time delays, we do not include it 
explicitly in the calculation. Instead, we will later consider 
its effects via an effective lens density profile slope, which 
will be the subject of detailed discussion. 

For the velocity function of early- type galaxies, w e 
adopt that derived from the SDSS data (IChoi et al.l ,20071) . 
which is fitted by a modified Schechter function of the fol- 
lowing form: 



dn 
dv 



v 

— ) exp 

v, 



P dv 
I>//3) V 



(4) 



Here, «*, a, /3)=(8.0 x 10- 3 /i 3 Mpc- 3 , 161 kms" 1 , 2.32, 
2.67). We do not include redshift evolution, i.e., we apply 
the mass distribution and velocity function to galaxies at 
any redshifts. Such a non-evolving model has been shown to 



reproduce th e abundance o f strong lenses in t he latest radio 
jChad ,20071 ) and optical dOguri et al l ,20Qg| ) lens 
although see, e.g.. lMitchell et al.l ( 20051 ) for an analysis that 



surveys, 



takes evolution into account. 



2.2 Lensing rates 



The lensing probability for a given source at redshift z s and 
with luminosity L is given by 



dO 



1 cdt 
dzi——(l 
dz\ 



Zl) 



dn dv 
'dvdB 



(Tie 



(5) 



u=u(0) 



where 6 is the image separation between multiple imagesQ z\ 
indicates the redshift of lensing objects, and dn/dv denotes 
the velocity function of lensing galaxies given by equation 
Q. Unless otherwise specified, we adopt the minimum im- 
age separation of 9 m in = 0.5" and maximum image sepa- 
ration of # max = 4". Although a non-negligible fraction of 
strong lenses have image separations larger than 4", such 
lensing is caused by groups or c lusters of g alaxies rather 
than a single massive galaxy (e.g- IOgurillioog ). In this work 
we focus on galaxy-scale lenses (for which simple isothermal 
density profiles are good models), for two reasons. Firstly, 
they do dominate the lens abundance, which we will need 
to be high in order to beat down the statistical errors on 
the cosmological parameters. Secondly, the image positions 
are easier to fit with simple models, reflecting the fact that 
massive galaxies are older, more relaxed systems than groups 



1 For strong lens systems withe multiple images more than two, 
we define the image separation by the maximum separation be- 
tween any pair of images. 



and clusters. On low mass scales, lenses with image separa- 
tion less than 0.5" were deemed to be too difficult to mea- 
sure, particularly in the ground-based surveys. Indeed, when 
we consider the yields of specific surveys below, we take into 
account the predicted image quality explicitly. 

The lensing cross section eri cns is determined by the 
structure of the lens potential. As discussed in i]2.1l we 
adopt an SIE plus external shear as the lens potential of 
each galaxy. To take account of th e magnification bias , we 
compute the biased cross section (|Huterer et all ,2005. ) as 
follows: 



trie 



I 



du d$/dL(L/fj,) 
~ d$/dL(L) ' 



(6) 



where d&/dL is the source luminosity function and /1 is 
the magnification factor. The integral is performed over 
the region where multiple images are generated. In prac- 
tice, we compute the biased cross section numerically using 
Monte Carlo sampling of multiple images. We solve the lens 
equation using the software package GLAFIC (M. Oguri, in 
preparation) . We compute the biased cross sections for dou- 
ble (two- images) , quadruple (quad; four-images), and naked 
cusp (three-images) lenses separately in order to study the 
image multiplicity. For double lenses, we place a condition 
that the flux ratio must be larger than 0.1 and perform the 
integral over the region that satisfies this condition. This 
is because asymmetric double lenses with large flux differ- 
ences are sometimes very difficult to locate in observations 
as a result of dynamic range problems caused by the brighter 
image. 

The choice of the magnification factor p depends on the 
way in which strong lens systems are identified in the sur- 
vey data. Bearing in min d the promisin g var iability selection 
techni ques suggested bv lPindorl (|2005l ) and lKochanek et alj 
(2006), we adopt the magnification factor of the fainter im- 
age for double lenses as /x. For quad and cusp lenses, we 
adopt the magnification factor of the third-brightest image, 
because the second-brightest image is sometimes located 
very close to the brightest image if it is produced near the 
fold catastrophe (where only one eigenvalue of magnifica- 
tion matrix is zero). The two images in a fold-configuration 
lens may be hard to deblend because of their small angular 
separation, which may result in failure to identify them as 
strongly lensed. Another reason is that the third brightest 
image tends to arrive first, a feature we will return to later 
(&■ 

Once the lensing probability for a given source is com- 
puted, it is straightforward to calculate the expected strong 
lens abundance. We obtain the differential number distribu- 
tion of strong lenses by integrating the lens probability over 
the source population: 

dN 

dz s 



d$ dV 
dM dMdz: P > 



(7) 



where p is given by equation © and M denotes the absolute 
magnitude of sources. The limiting absolute magnitude is 
simply converted from the magnitude limit of the survey 
considered. The volume factor dV/dz s is given by 



, cdt 



QDt — (l + z s y 



(8) 



dV_ 

dz s E dz s 

with being the solid angle corresponding to the survey 
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Finally, the total number of strong lenses is 

dN 



N = dz 



' dz K 



(9) 



Note that we did not impose any restrictions on the lu- 
minosity of the lensing galaxies. If a lens galaxy is much 
brighter than lensed images, such lens systems may prove 
difficult to identify from the imaging data alone. This ef- 
fect becomes particularly important for lensing o f faint 
sources (e.g., lKochaneklll996l ; IWvithe fc Loebll2002h . How- 
ever, since we are considering cadenced surveys where the 
strong lens search will make use of differenced survey images 
l|Kochanek et alj|2006h . we assume that those systems with 
bright lensing galaxies can still be successfully identified be- 
cause the lens galaxy component (as well as any lensed host 
galaxy light) should be cleanly subtracted. In practice, a 
bright lens galaxy will act as an additional source of noise; 
we neglect this here, since many of the surveys we consider 
are ground-based and will have background-limited images. 




-26 -28 
M [mag] 



2.3 Number density of quasars 

We adopt the standard double power-law for the quasar lu- 
minosity function. Specifically, we adopt the following para- 
metric form: 



QSO 



4>, 



dM 



lQ0A(a+l)(M-M.) _|_ lQ0A(fi+\)[M-M,) ' 



(10) 



Here M refers to the absolute i-band magnitude of quasars. 
We fix the faint end slope of f3 — —1.45 and the bright end 
slope of a = —3.31 which was obtained in the comb ined 
analysis of the SDSS and 2dF (jRichards et al.ll2005h . At 
z > 3, however, we modify the bright end slope and adopt 
the shallower slope of a = — 2.58 which was suggested by 
observations l|Fan et al.ll200lh . 

In addition to this change of slope, we model the red- 
shift evolution of the luminosity function as pure luminosity 
evolution. Specifically, the break absolute magnitude M* is 
described as 



-20.90 + 5 log ft- 2.5 log f(z) 



/(*) 



s Cz (l + e' 



(11) 



(12) 



where the zero-point value of M* , as well as the normalisa- 
tion of the lu minosity functio n cfr* = 5.34 x 10~ 6 /i 3 Mpc -3 , 
is taken from [Richards et al.l l)2005h . The par ametric form 
of the redshift evolution, f(z), is taken from iMadau et all 
l| 19991 ) with some modifications to achieve a better match to 
the observed evolution of the quasar luminosity function. We 
then determine the values of the parameters £, £, and z* to 
reproduce the quasar luminosity function fro m z = to 5 ob- 
taine d from the SDSS Data Release 3 (DR3, I Richards et al.l 
l2006t ). We find the best fit values to be ((, f , z,)=(2.98, 4.05, 
1.60), which we adopt throughout the paper. In Figurefl] we 
compa re our best-fit lumino sity function with the SDSS re- 
sult of iRichards et al.l (120061 ). 

In order to convert the i-band absolute magnitude 
to the apparent magnitude, we need the appropriate K- 
correction. We adopt that presented by IRichards et al.l 
<|2006l ). 



Figure 1. Our model of the i-band quasar luminosity function 
(cq. 1101 ). compared wi th the observed quasa r luminosity function 
from the SDSS DR3 llRichards et alj|2006l) . We present results 
at three different redshifts, z = 1.25, 2.01, and 3.75, in order to 
show that our model successfully reproduces the abundance of the 
SDSS quasars for a wide range of redshifts. Luminosity functions 
for z = 1.25 (z = 3.75) are shifted by —1 dex (+1 dex) to avoid 
overlap. 



2.4 Number density of supernovae 

The rate of occurrence of SNe is closely related to the star 
formation rate. In this p aper we adopt the cosmic s tar for- 
mation rate presented bv lHopkins fc Beacoml (|2006l ). 



, . (0.0118 + 0.08z)h. llr _ u , _ 
PSfr(*)= 1 + (3/3.3)5.3 M P C 



(13) 



which is essentially the best-fit to the observed star 
formation rates assuming the initial mass function of 
iBaldrv fc Glazebrookl (|2003h . 

Recent studies of the host galaxies of SNe la suggest 
that delay times (to) of SNe la cannot be described by a 
single value. It has sometimes been assumed that SNe la are 
drawn from two populations, "prompt" (tracing star forma- 
tion r ates) and "old" (tracing stellar masses of host galaxies) 
SNe JScannapieco fc Bildstenl 120051; iMannucci et al.l 120051 . 



Sullivan et al.ll2006bl ; lAubourg et al.l l2008). However, 



l200d 

such a two-component model may fail to ex plain the ob- 
serve d low high-redshift SNe la rates (e.g., iDahlen et al.l 
2008). It is naturally expected that more realistic mod- 
els have widely distributed trj, because theoretically the 
delay time should depend on various parameters such as 
the mass of the companion star and the metellicity (e.g., 
iKobavashi fc Nomotol [20091 . and reference s there in). In this 
paper, we use the result of iTotani et al.l |2008) who con- 
strained the distribution of to assuming a power-law: 



/(to) oc *d (to > O.lGyr). 



(14) 



Since they found no strong evidence for the presence of the 
"prompt" (to < O.lGyr) component, we do not consider it. 
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Figure 2. The evolution of SN rates for type-la (lower line) 
and core-collapse (i.e., the sum of Ib/c, IIP, IIL, and Iln; up- 
per line) SNe adopted in this paper. See text for details of the 
model. Filled s quare s are recent m e asurem e nts of SN la ra t es from 
Hardin et ail ||200(]|1. iPain et all J2002T). iDahlen et all d2004) , 
Blanc et aD <2004h. iNeill et al.l i2006h. IPoznanski et ail 120071). 
Kuznetsova et ail j20Qg|'). IBotticella et al.l [120081), iDildav et al] 
l|2008h . iHoresh et alf " il200Sft . and IDahlen et al.l l|200Sl ). whereas 
open squar es are recent meas ur ements of core - collap se SN 
rates from IDahlen et ail J2004li IBotticella et ail d2008Tl . and 
Bazin et al. I [|2009ri . Errors indicate lcr statistical errors, and do 
not include any systematic errors. 



We can then compute the SN la rate as 



nia(2) = JJCl 



/ 1 PSFR[z(t - t D )]f(t D )dtz 



(15) 



The factor Ci a = 0.032A/q can be computed from the 
stellar mass range of 3Mq < M < 8Ma for SNe la an d 
the initial mass function of iBaldrv fc Glazebrookl l|2003h . 
We assume a canonical efficiency of 77 = 0.04 (see also 
IStrigari et aill2005l : [Hopkins fc Beacoml 120061 ). 

On the other hand, core-collapse SNe are thought to be 
related more directly to the star formation rate. They are 
associated with the death of massive stars whose life time 
is significantly shorter than the typical cosmological time 
scale, which suggests that the SN rate of core-collapse SNe 
is simply proportional to the cosmic star formation rate: 



ricc(z) = C cc psFB.(z), 



(16) 



where C cc = 0.0132A/Q 1 was derived assuming that core- 
collapse SNe are produced in the stellar mass range 8Mq < 
M < 50M(7), together with the initial mass function 
of IBaldrv fe Glazebrookl l|2003h . Following lOda fe Totaml 
(2005), we adopt constant relative proportions for four sub- 
classes of core-collapse SNe. Then SN rates are written as 



nx(z) = fxn cc (z), 



(17) 



wher e ,fi b / c = 0-23, /hp = 0- 30, /iil = 0.30, and /n n = 
0.02 IDahlen fc Franssonl \l99^ ) . Note that we do not con- 
sider SN1987A-like SNe that are very faint compared with 



the other core-collapse SNe. Recent analysis suggests that 
the fraction of type -HP SNe may actually be much higher 
l|Smartt et alj|2009l ). 

In Figure [2] we compare our model of SN rates with 
recent measurements. We find that our model is reasonably 
consistent with the observed SN rates, particularly given the 
additional systematic errors they involve. The Figure indi- 
cates that our model of SN rates is uncertain at the factor of 
~ 2 level, suggesting that our predicted numbers of lensed 
SNe will be similarly uncertain. In particular, our model 
appears to slightly over-predict the core-collapse SN rate, 
though it is possible that current observations may miss a 
fraction of core-collapse SNe given their wide range of in- 
trinsic luminosities. The uncertainty should be even larger 
at redshift z > 2 for type-la SNe, and z > 1 for core-collapse 
SNe, where no measurements of SN rates have been ob- 
tained. 

We need not only SN rates, but also the brightness dis- 
tributions of SNe in order to make predictions of the lensed 
SNe abundance. In what follows the magnitude of SNe refers 
to the peak magnitude, i. e., the magnitude wh en the SN is 
brightest. Again following [Oda fc Totanll |2005l '), wc assume 
the abso lute magnitudes of SNe a re Gaussian-distributed 
(see also lYasuda fc Fukugital l2010t ). With this assumption 
the luminosity function (in terms of B-band absolute mag- 
nitude M) can be written as 



dM 



nx(z) 



(1 + z) v^crx 



exp 



(M - Mi) 2 



2al 



(18) 



where (Mf a , M : * b/c , Mf IP , M„ L , Mf Ia )=(- 19.06, -17.64, 
-16.60, -17.63, -18.75) (for h = 0.72) an d (oia, <7i b /c, ^iip, 
am, criin)=(0.56, 1.39, 1.12, 0.90, 0.92) (jRichardson et ail 
l2002f ). Note that the luminosity function of SNe (eq. [18]) 
differs from that of QSOs (eq. [10]) in that the former is in 
fact the number rate (number per unit time). Thus the fac- 
tor of (1 + z) -1 is introduced to account for the cosmological 
time dilation. 

We convert B-band absolute magnitudes to apparent 
magnitudes in the i-band by computing K-corrections using 
various SN templa te spectra. W e adop t the sp ectra at the 
peak presented by iNugent et ail (l2002f) fo r la, [l evan et ail 
(|2005h for Ib/c, and iGilliland et ail (|l999T ) for IIP, IIL, and 
Iln. 



3 STRONG LENSES IN VARIOUS SURVEYS 

In this section, we predict the number of strongly lensed 
QSOs and SNe for a selection of ongoing and planned sur- 
veys, using the model described in detail in 32] 



3.1 The expected number of lenses as a function 
of survey depth 

First, we explore how the number of lenses detected increases 
with survey depth. In Figure [3] we plot the number of lensed 
QSOs (in a half-sky survey) as a function of i-band limiting 
magnitude % m . The slope of these number counts is fairly 
shallow, particularly at iij m > 21, which suggests that the 
survey area is much more important than the survey depth 
when trying to discover many strongly lensed QSOs. The 
lensing rate is ~ 10~ 3,5 , and does not depend very much on 
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Figure 3. The expected number of lensed QSOs as a function 
of the i-band limiting magnitude iu m . A fiducial survey area of 
f! = 20, 000 deg 2 is assumed. The number of non-lensed QSOs 
is also shown for reference. The lower panel shows the ratio of 
the number of quad (four-images) or naked cusp (three-images) 
lenses to the total number of lenses, as a function of i\i m . 



in m due to the conflicting effects of increasing mean QSO 
redshift and decreasing magnification bias. We note that the 
lensing rate is lower than observe d in the Cosmic Len s All- 
Sky Survey (CLASS), ~ 10" 2 ' 8 dBrowne et alj|2003h . pre- 
sumably because of the quite different magnification bias it 
involves (for instance, in the CLASS the total magnification 
factor is used for the magnification bias). The recent optical 
lens survey, the SDSS Quasar Lens Search (SQLS), has ob- 
tained a lensing rate of ~ 10~ 3 3 , which is more consistent 
with the calculation above (see llnada et al. 2008). The frac- 
tion of quad lenses decreases from ~ 30% for in m = 18 to 
~ 10% for iii m = 2 8, which is roughly consi s tent with pre- 
vious calculations l|Rusin fc Tegmarkl l200ll ; iHuterer et al.l 
l2005l ; IOgurill2007bl ). The small fraction (~ 10" 3 ) of naked 
cusp lenses indicates that only very wide-field surveys will 
be able to locate such rare image configurations. We again 
note that our calculation is applicable only to galaxy scale 
lenses; naked cusp lenses are much more common at clus- 
ter scales, wh ere the radial density profiles of the lenses 
are shallower (|Qguri fc Keetonl |2004 iMinor fc KaplinghatJ 
2008). We also note that our model lens galaxies are all 
spheroids, with correspondingly low ellipticity. The lower 
mass, disky lenses, that will make up a small minority 
(< 20%, see £|2.1[) of any survey's yield, may be expected to 
be more elongated and hence prese nt a higher naked cusp 
fraction l|Keeton fc Kochanekll 19981 ). 

In Figure|4]we plot the expected numbers of lensed SNc. 
The slope of the number counts is steeper than it was in the 
case of the lensed QSOs, suggesting that the survey depth is 
more important for lensed SNe than it is for lensed quasars. 
Lensed core-collapse SNe are more abundant than lensed la 
by a factor of several. We find that the fraction of quad 
systems changes very rapidly as a function of the limiting 



Figure 4. The expected number of lensed SNe as a function of the 
i-band limiting peak magnitude, i pe ak,llmi f° r one year monitoring 
of a fiducial f2 = 20,000 deg 2 survey area. Upper curves indicate 
numbers of normal unlensed SNe, for reference. Solid lines show 
expected numbers of SNe la, whereas dashed lines show the num- 
bers of core-collapse SNe. The thin lines at > 24 mag show the 
expected numbers of SNe as a function of ff pe ak lim rather than 
*peak,limi as relevant for a near-infrared optimised space-based SN 
survey. The lower panel shows the fractions of quad and naked 
cusp lenses. 



magnitude. For shallow SN surveys, with, say, limiting peak 
magnitude i pea k,lim J$ 22, the strong lens sample would be 
approximately 50% quads. The high quad fraction is due to 
very large magnification bias inherent to such bright lenses. 
For deeper surveys of i pea k,iim > 26, however, the fraction of 
quad lenses takes a typical value of 10 — 20%. The fraction 
of cusp lenses remains low, ~ 10 -3 ' 5 . 

The number of (lensed) SNe does not increase very 
much beyond i pca k.iim ~ 24. This is mainly because of the in- 
creasing values of the K-correction at high-redshifts. Since 
SN spectra have UV cutoff at ~ 3000 A, the i-band K- 
correction rises quickly at 2 > 1.7, which significantly de- 
grades the detectability of such high-redshift SNe in i-band 
images (this effect has long been a source of frustration for 
ground-based supernova cosmology projects). To illustrate 
this point, in Figure [4] we also plot the number of SNe as a 
function of //-band limiting peak magnitude, i/ pea k,iim- It is 
clear that the expected number of lensed SNe is much larger 
in _ff-band than i-band for deep monitoring. The steeper 
number counts in //-band indicates that the effect of K- 
correction is much less important, illustrating the impor- 
tance of near-infrared imaging for detecting high-redshift 
SNe (lensed or otherwise). 



3.2 The expected number of lenses in individual 

surveys 

Next we consider strongly lensed QSOs and SNe in several 
specific time-domain surveys. Each survey is characterised 
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Table 1. Properties of various time-domain surveys. These are fiducial numbers based on assumptions we describe in the text; given are 
survey area Q, 10a point source limiting magnitude irim for one visit, one year, and for the final stacked survey image, corresponding 
assumed median image quality #pgp (seeing FWHM, in units of arcseconds), cadence (in days), season length t BeaBon (in months), survey 
length t surV cy (in years) and the effective survey duration t e ff = i survcy • t seaaon /12 (also in years). 



Survey 


n 


kirn (visit) 


«lim(year) 


him (total) 


0PSF 


cadence 


^season 


^survey 


*cff 




[deg 2 ] 








[arcsec] 


[days] 


[months] 


[years] 


[years] 


SDSS-II 


300 


21.3 


22.9 


23.5 


1.4 


5 


3 


3 


0.75 


SNLS 


4 


24.3 


26.3 


27.1 


0.9 


4 


5 


5 


2.1 


PS1/37T 


30000 


21.4 


22.2 


22.7 


1.0 


30 


2 


3 


0.5 


PS1/MDS 


70 


23.3 


25.2 


25.8 


1.0 


4 


6 


3 


1.5 


DES/widc 


5000 


23.6 


24.0 


25.1 


0.9 


30 


1 


5 


0.4 


DES/deep 


6 


24.6 


26.1 


27.0 


0.9 


5 


3 


5 


1.3 


HSC/wide 


1500 


24.9 


24.9 


25.5 


0.75 







3 


0.0 


HSC/deep 


30 


25.3 


26.6 


26.6 


0.75 


5 


2 


1 


0.2 


JDEM/SNAP 


15 


27.1 


. . b 


29.7 


0.14 


4 


12 


1.3 


1.3 


LSST 


20000 


23.3 


24.9 


26.2 


0.75 


5 


3 


10 


2.5 



a The HSC/wide survey is zero cadence: with only one observing epoch per year it would not be suitable for useful lensed SN detection. 
'The JDEM/SNAP survey is assumed to be undertaken as a single 1.3 year observing campaign, so does not offer the yearly stack 
detection strategy. In practice lensed quasars would be detectable simply from the resolved image geometries in any case. 

Table 2. The expected numbers of lensed QSOs in various time-domain surveys. The numbers are either for detections only, using 
stacked images within each survey year, or for well-measured time delays using each visit's image. Some wide-field surveys do not have 
high enough cadence to measure time delays, and therefore those surveys are useful only for detecting lensed QSOs (see text for more 
details). We adopt the minimum image separation (9 m ; n = (2/3)$pgp for all surveys. The numbers of non-lensed QSOs detectable in the 
surveys are also shown for reference. Percentages in parentheses indicate the fraction of quad lenses. 



QSO (detected) QSO (measured) 



Survey 


N i 

1 " nonlcns 


-^lcns 


N i 

1 ' nonlcns 




SDSS-II 


1.18 


X 


10 5 


26.3 (15%) 


3.82 


X 


10 4 


7.6 (18%) 


SNLS 


9.23 


X 


10 3 


3.2 (12%) 


3.45 


X 


10 3 


1.1 (13%) 


PS1/37T 


7.52 


X 


10 6 


1963 (16%) 










PS1/MDS 


9.55 


X 


10 4 


30.3 (13%) 


3.49 


X 


If) 4 


9.9 (14%) 


DES/wide 


3.68 


X 


10 6 


1146 (14%) 










DES/deep 


1.26 


X 


10 4 


4.4 (12%) 


6.05 


X 


10 3 


2.0 (13%) 


HSC/wide 


1.76 


X 


10 6 


614 (13%) 










HSC/deep 


7.96 


X 


10 4 


29.7 (12%) 


4.30 


X 


10 4 


15.3 (13%) 


JDEM/SNAP 


5.00 


X 


10 4 


21.8 (12%) 


5.00 


X 


10 4 


21.8 (12%) 


LSST 


2.35 


X 


10' 


8191 (13%) 


9.97 


X 


10 6 


3150 (14%) 



Table 3. The expected number of detected lensed SNe (type la and core collapse) in various time-domain surveys. We adopt the 
minimum image separation 6 mln = (2/3)0psf for all surveys. The numbers of non-lensed sources detectable in the surveys are also shown 
for reference. Percentages in parentheses indicate the fraction of quad lenses. For lensed SNe, we adopt the peak magnitude limit of 
*peak,lim = ilim — 0.7 in actual calculations so that the lightcurves of lensed SN images can well be traced. 



Survey 




SN (la) 


SN ( 


cc) 


Note 


■^nonlens 


-^lens 


-^rionlens 


-^lcns 


SDSS-II 


4.34 x 


10 2 


0.003 (54%) 


1.09 x 10 3 


0.01 (40%) 




SNLS 


7.52 x 


10 2 


0.03 (24%) 


1.44 x 10 3 


0.05 (26%) 




PS1/37T 


3.34 x 


10 4 


0.28 (53%) 


8.23 x 10 4 


0.97 (39%) 


detections only 


PS1/MDS 


2.93 x 


10 :! 


0.09 (32%) 


6.05 x 10 3 


0.16 (30%) 




DES/wide 


8.30 x 


10 4 


2.7 (29%) 


1.62 x 10 5 


4.9 (29%) 


detections only 


DES/deep 


8.95 x 


10- 


0.04 (22%) 


1.80 x 10 3 


0.07 (24%) 




HSC/deep 


1.10 x 


10 :i 


0.06 (18%) 


2.56 x 10 3 


0.13 (21%) 




JDEM/SNAP" 


1.36 x 


10 4 


2.9 (13%) 


5.39 x 10 4 


12.0 (18%) 




LSST 


1.39 x 


10 6 


45.7 (32%) 


2.88 x 10 6 


83.9 (30%) 





a Instead of the i-band, we adopt an .ff-band magnitude limit of Hy lIn = 26.8 to predicted the number of (lensed) SNe, since in practice 
the detection in space will be done in the near-infrared to optimise the number of high-redshift SN sources. 
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by the survey area (fl) and AB magnitude limit iu m , de- 
fined by the i-band 10a detection limit for point sources. 
In some cases, we estimated an approximate lOcr limit from 
the 5cr detection limit quoted in the reference by subtracting 
0.7 mag. In some cases we also made rough conversions be- 
tween image qualities, in order to represent accurately the 
conditions currently expected to be achievable. These ap- 
proximations were confirmed to be accurate at the 5% level 
using various exposure time calculation toolsjf] For lensed 
SNe the interpretation of the magnitude limit is not obvi- 
ous because the brightnesses of SNe change drastically with 
time. We set the condition that (lensed) SNe will be detected 
if the peak magnitude is more than 0.7 mag brighter than 
the magnitude limit of each visit. This requirement is chosen 
to ensure that the lightcurve of each SN is traced reasonably 
well, and is con sistent with the 25a peak brightness detec- 
tion required by I Young et alj (|2008T ). In addition, for lensed 
SNe we need to specify the effective survey duration t e g to 
predict the total number of lensed SNe. Depending on the 
survey duration and design, a fraction of lensed SN events 
may be missed simply because one of the multiple images 
appears outside the survey tim e; in this paper w e do not 
consider this "time delay bias" (|Oguri et al.l 120030 for sim- 
plicity. We also adopt different minimum image separation 
#min for different surveys, depending on their typical seeing 
sizes. Specifically we assume f? m i n = (2/3)#psf, where #psf 
is a typical seeing FWHM of each survey, which is empiri- 
cally consistent w ith with previous optical lens surveys (e.g., 
Ilnada et alj|200Sf ). 

Finally, we note that detecting a time-variable lens, and 
being able to measure its time delay, are two very different 
things. Since q uasars are marked out by their long term 
variability (e.g., iKellv et al, I l2009l . Schmidt et al, in prepa- 
ration), a good strategy for detecting lensed quasars may 
be to stack images within each survey year, and then look 
for variability between years. For this reason we consider 
the limiting magnitude of a survey per year's observation 
for detecting quasar lenses. Measuring time delays requires 
dense sampling: we assume that a cadence of < 10 days and 
a season length of ^ 3 months are the minimum needed to 
allow time delay estimation, based on results of pre vious lens 
monitoring campaigns fe.g. iFassnacht et aJj|200"3 ). We note 
that quantifying time delay accuracy as a function of pho- 
tometric precision, season length, and cadence would make 
a very valuable future study. Detecting useful lensed SNe 
will also require high cadence survey data, since for these 
lenses there is no option of building up time delay precision 
in subsequent seasons. 

We now briefly describe ongoing and future time- 
domain surveys which we consider in this paper. Survey pa- 
rameters for these surveys, including survey areas, limiting 
magnitudes, and survey durations, are listed in Table [1] 

• Sloan Digital Sky Survey-II Supernova Survey (SDSS-II 
SN). This survey was designed to locate relatively low- 
redshift (z ~ 0.05 — 0.35) supernovae in the so-called "stripe 

2 Specifically, we used the CFHT MegaCam facility "DIET," 
http : //rpm. cf ht . hawaii . edu/~megacam/diet/DIET . rpm, 
and also the LSST ETC at 

http : //dls .physics .ucdavis . edu: 8080/etc4_3work/ 
servlets/LsstEtc .html 



82" region of the SDSS (|Frieman et al.ll2008f ). It involved 
monitoring an area of 300 deg 2 with five broad-band filters; 
the limiting magnitude at each epoch is in m = 21.3. The 
survey consisted of three three-month observing seasons, 
corresponding to an effective survey duration for the SN 
search of t e g = 0.75 yr. We assume a typical seeing of l'.'4, 
and a monitoring cadence of 5 days, giving 18 visits per 
field per year. 

• Supernova Legacy Survey (SNLS). Running for 5 years, 
with roughly five-month observing seasons for a total effec- 
tive survey duration of 2.1 years, the SNLS is the cadenced 
part of the CFHT legacy survey, designed for discovery, 
classification and monitoring of supernovae at intermediate 
redshift. The much larger "Wide" part of the survey has 
only static images that do not (straightforwardly) enable 
a time delay lens search. The "Deep" fields cover 4 deg 2 . 
In calculating the i-band magnitude limit ir im = 24.3, we 
assumed (consta nt) median seeing / . / 9 , and 1800s exposure 
time per epoch (jSullivan et al . 2006a). We also assume a 

4- day cadence, leading to 38 visits per field per year. 

• Pan-Starrs 1 (PS1). The first 1.8-m telescope of the 
Pan-STARRS project is carrying out two surveys, the 
Medium Deep Survey (MDS) and the 37r survey. Both are 
cadenced, but at different rates. The MDS covers 70 deg 2 , 
and we assume that each of its fields will be observed 
for three six-month seasons with high (typically 4-day) 
cadence. We anticipate the 3ir survey covering 30,000 deg 2 
at lower (typically one month) cad ence, and assume t hree 
three-month seasons (we follow lYoung et al. I 120081 , m 
anticipating using all the filters in the SN identification). 
We conservatively assume a median image quality of l'.'O, 
and use this to estimate 10cr i-band single-visit limiting 
magnitudes for each survey of 23.3 (MDS) and 21.4 
(37r), extrapolating approxim ately from the figures give n 
in the design reference mission (jChambers fc DennaulboOSf ). 

• Dark Energy Survey (DES). To be carried out with 
a new wide-field camera installed on the Blanco 4-m 
telescope at CTIO, this survey, to start in 2012, is also 
planned to have two components, a wide a nd a deep 
(|The Dark Energy Survey Collaboration! l2005h . As with 
PS1, only the deep fields will have sufficient cadence for 
time delay estimation, but the wide fields will again allow 
detection of time-delay lenses by their variability. While 
the details of the surveys are still being decided, we adopt 
the following fiducial estimates: a deep survey covering 
6 deg 2 for five years (in nearly three-month seasons, with 

5- day cadence, giving 16 visits per year), and a wide survey 
covering 5000 deg 2 for fi ve years, with two observing epochs 
per year (lFriemani r2009). We take the effective season length 
of the wide survey to be one month, in order to estimate 
the number of lensed supernovae seen: this is about the 
period over which 2 images might be simultaneously visible. 
We assume iu m — 24.6 per epoch in the deep fields, and 
in m = 23.6 per epoch in the wide fields (J. Annis, priv. 
comm.). 

• Hyper Suprime-Cam (HSC). The HSC is a next- 
generation prime-focus camera for the Subaru 8.2-meter 
telescope. The field-of-view is 1.5 degree in diameter, and 
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the expected image quality is about 0'.'75 across the entire 
field-of-view. The survey plan with the HSC is not yet 
finalised and still very uncertain, but one of the three survey 
components ("deep") has both area and depth enough to 
locate many high-redshift supernovae. Here we tentatively 
assume its survey area to be 30 deg 2 and the limiting 
magnitude to be iu m = 25.3, with a single two month 
season (which therefore may not be sufficient to measure 
time delays between multiple images). At an assumed 
cadence of 5 days, each field would be observed 12 times. 
Another component (the "wide") may cover 1500 deg 2 , with 
the imaging again being built up one exposure per year, 
for two or three years. As with PS1/37T and DES/wide, we 
flag this survey as useful for detection but not time delay 
measurement. We do not consider the "ultra-deep" survey 
here as its survey area will be rather small. 

• Joint Dark Energy Mission (JDEM). The design of 
this space-based observatory is still very uncertain, but 
we can assume that if such a facility is built with high 
redshift supernovae in mind, it will have an infrared- 
optimised imager. In this work we compute, as an example, 
predictions for a JDEM that has indeed been optimised 
for supernova work, and assume the specifications of the 
Supernova / Acceleration Probe (SNAP) as our fiducial 
telescope l|Aldering et alj [20041 ): we assume a 15 deg 2 sky 
area observed to a depth of i !im = 27.1 and Hu m = 26.8 
per epoch at high cadence (4 days) for 1.3 years (where 
we have again corrected to get 10a limits). Note that we 
only consider in this work the SN survey of this strawman 
JDEM; a weak lensing survey in space should find large 
numbers of lensed quasars even if the imaging is not 
cadenced at all, as the single epoch high resolution imaging 
alone wil l allow identificatio n of multiple image systems 
(see, e.g.. lKuhlen et al.ll2004 for some estimates). However, 
time delay measurements would be possible only in the SN 
fields. 

• Large Synoptic Survey Telescope (LSST). This 6.7-m 
effective aperture dedicated survey telescope is designed 
to survey 20,000 deg 2 of sky for ten year s . The "universal 
cadence" survey outlined by llvezic et al.l l|2008h sees each 
field being observed in 6 filters at typically 5-day cadence 
in the r and i-bands. A typical field's observing season 
would last three months or more (with the redder filters 
sampling a longer season), giving some 20 exposures per 
field per filter per year. We estimate that the 10<r limiting 
magnitude in CI? 75 seeing will be aroun d 23.3, and adopt 
this as our fiducial single-epoch depth, llvezic et al.l note 
that "mini surveys" could be carried out with LSST on 
smaller sky patches to greater depth; we do not consider 
any such survey here, but note that it would provide higher 
accuracy time delays on an inevitably smaller number of 
lenses. 



Tables [2] and summarises the expected numbers of 
lensed QSOs and SNe in each survey. We find that the LSST 
will find the largest number of lenses with measureable time 
delays, for both quasars and SNe. The wide surveys of the 
HSC, DES and PS1 projects will also discover large num- 
bers of lensed quasars thanks to their wide field-of-view, 



but they will not have sufficient time domain sampling to 
measure the time delays between quasar images. The cor- 
responding deep surveys, or a spaced-based JDEM/SNAP 
survey, will provide time delay measurements, but only for 
~ 60 lensed quasars, because of their small combined field- 
of-view. The LSST is unique because it has both the wide 
field-of-view to cover many quasars and the frequent time 
domain sampling for monitoring of lens systems, allowing 
us to obtain time delay measurements for ~ 3000 lensed 
quasars. It will be quite challenging to find lensed SNe in 
ground-based surveys before LSST, as the expected number 
of lensed SNe discovered is only ~ 1. In contrast, the very 
deep monitoring by a JDEM/SNAP survey would allow us 
to locate ~ 15 lensed SNe in total, despite the relatively 
small field of view of 15 deg 2 . The LSST will also be power- 
ful for detecting lensed SNe; our calculation indicates that 
the LSST will discover ~ 130 lensed supernovae during its 
10-year survey. 



4 GENERATING MOCK CATALOGUES: A 
MONTE-CARLO APPROACH 

Here we compute the expected number of strong lenses, and 
their distributions with lens and source property, in a differ- 
ent way from the standard technique described in ffj?] Essen- 
tially it is a semi-analytic technique based on Monte-Carlo 
realisations of lens and so urce populations, and is an exten- 
sion of the method used in lOguri fc Blandfordl l|2009f) to pre- 
dict the all-sky distribution of large Einstein radii produced 
by massive clusters. In essence, we perform the necessary 
integrals by Monte-Carlo integration, and keep the sample 
points drawn from the target distribution. The goal is to 
produce mock catalogues of lenses for a given survey, which 
can then be used to investigate the feasibility of various sci- 
ence projects. 

4.1 Technique description 

The specific procedure is as follows. First we randomly gen- 
erate a list of sources in a given survey area to a given 
depth, according to the adopted source luminosity function. 
To allow for the magnification by strong lensing, the limit- 
ing magnitude of the source catalogue needs to be deeper 
than the actual magnitude limit of the survey. We adopt 
3 — 3.5 mag deeper limit than the survey magnitude limit 
for generating the source catalogue. This is sufficient, partic- 
ularly given our use of the fainter image magnification factor 
in computing the magnification bias. Each source realisation 
contains information on the redshift and apparent magni- 
tude (without lensing magnification) only. Since we ignore 
any spatial correlations between sources, the whole source 
population is characterised just by a surface number density. 
We then randomly generate a list of lenses according to the 
model described in §2.1\ For each mock lens galaxy, we con- 
sider the rectangular region with an area of [8#Ein(max)] 2 , 
where 6 l Ei n (max) is given by equation with D\ s /D s = 1, 
randomly distribute source objects using the pre-generated 
catalogue, and for each source solve the lens equation to see 
if it is multiply-imaged. Again, the lens equation is solved 
using the software glafic (M. Oguri, in preparation). For 
multiply-imaged sources, we check whether they satisfy the 
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Figure 5. The predicted redshift distributions of lensed QSOs 
(upper panel) and SNe (lower panel). The histograms are dis- 
tributions calculated from the mock lens catalogue constructed 
using the Monte-Carlo technique, whereas the solid curves are 
analytic results from equation 



detection criteria (the survey limiting magnitude, image sep- 
aration, and flux ratio), and record only lens systems which 
survive these criteria. In this way we can generate mock 
catalogues of lensed sources for a particular survey. 

Although the new Monte-Carlo technique tends to 
be more time-consuming than the traditional analytic ap- 
proach, it has several important advantages. (1) In this ap- 
proach we can add more parameters to the lens model with 
little additional computational cost. In the traditional ana- 
lytic approach, adding additional lens parameters requires 
an additional integral, particularly if the distributions of the 
parameters are not independent of each other. This exponen- 
tial growth in run time with model complexity is in marked 
contrast with the Monte-Carlo approach, in which the CPU 
time required scales approximately linearly with the number 
of model parameters. (2) It is straightforward to include even 
quite complicated lens selection functions, e.g., those depen- 



Figure 6. Same as Figure \E\ but the redshift distributions of 
lensing galaxies are shown. 



dent on many observables such as flux ratios, image separa- 
tions, time delays, and the properties of lensing galaxies. (3) 
From the resulting mock catalogue we can examine the ex- 
pected distributions of lens and source parameters quite eas- 
ily. This is important because the population of strong lenses 
generally differs from the population of general lens and 



source populations (e. g., Oguri et al. 2005): Moller et all 



120071 : iRozo et all [2007: Mandclba uin et alJlioolih . and yet 



we would very much like to use the former to understand 
the latter. 



4.2 Mock lens catalogue 

As a specific example, in this paper we present the mock 
catalogue of lensed QSOs and SNe expected for the baseline 
survey planned with LSST (see Table [1] for the survey pa- 
rameters). In practice, the catalogue is 5 times (QSOs) or 
10 times (SNe) over-sampled in order to reduce shot noise. 
The catalogue contains all the necessary information for each 
mock lens, including the properties of lenses and sources, 
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Figure 7. Upper panel: The probability distribution of the lens 
galaxy ellipticity e for double (solid) and quad (dotted) lenses, 
obtained from the mock lens catalogue. Lower panel: The prob- 
ability distribution of 8 e — 7 , where 8 e is the position angle of 
the lensing galaxy and 7 is the position angle of external shear. 
Since we assumed random distributions for 9 e and $ 7 , the original 
unbiased (non-lens) population has a flat distribution of 9 e —8-y. 

image positions, magnifications, and time delays between 
image pairs 

First, we examine the redshift distributions, to check 
the validity of the Monte-Carlo technique. Figure[5]compares 
the source redshift distribution of the mock lenses with the 
analytic result obtained from equation ©. The reasonable 
agreements of the distributions assure us that the Monte- 
Carlo technique is indeed feasible and reliable. The total 
number of strong lenses in the mock catalogue is 3132 ± 25 
for lensed QSOs and 122 ±4 for lensed SNe (the errors refer 
to the lcr Poisson error), which is again in good agreement 

3 The mock lens catalogues described here are available at 
http : //kipac-prod. Stanford. edu/collab/research/lensing/mO' 



with the analytic result shown in Tables [2] and [3] In Figure [6] 
we show lens redshift distributions which are again in good 
agreement. 

As discussed above, the mock lens catalogue is useful in 
studying the lensing bias, i.e., the difference of the lens pop- 
ulation from the normal population. As practical examples, 
in Figure [7] we show the distributions of the lens ellipticity, 
and of the alignment between the lens galaxy and external 
shear, for double and quad lenses separately. We find that 
the lensing galaxies of quad image systems are Ae ~ 0.1 
more elliptical than th ose of double len ses. This result is 
broadly consistent with lRozo et al.l l|2007l ) who showed such 
a difference using the analytic approach. It appears that 
there is no significant difference in the distributions of ex- 
ternal shear between double and quad lenses. However, we 
find that the relative orientation of lensing galaxies and ex- 
ternal shear is biased: for quad lenses the direction of exter- 
nal shear is more likely to be aligned with the major axis of 
the lensing galaxy, whereas for double lenses the direction 
of external shear is more likely to be aligned with the minor 
axis. This bias can easily be understood by the well-known 
approximate degenera cy between galaxy ellipticity and ex- 
ternal shear (see, e.g. jKeeton et al!ll997l ). 



4.3 Distributions of time delays 

The predicted distributions of time delays are important in 
order to assess the measurability of time delays from fu- 
ture time-domain surveys. In Figure[8]we show the expected 
distributions of time delays and magnitude differences for 
strongly lensed image pairs, derived from our LSST mock 
lens catalogue. For double lenses, typical time delays are 
~ 1 — 3 months. Lensed SNe have somewhat shorter time 
delays on average than lensed QSOs, partly because of their 
lower redshifts and the stronger effect of the magnification 
bias. In most cases the brighter images arrive first, as is ex- 
pected from simple spherical mass models, but in some cases 
the arrival order is inverted, i.e., the fainter images arrive 
first. In either case the detection of first arrival images is as- 
sured because of our use of the fainter image for computing 
the magnification bias. 

For quad lenses the situation is more complicated be- 
cause there are three independent image pairs. However, in 
Figure [H] we can clearly see a general trend that the first 
arrival image tends to be fainter than the second and third 
arrival images (Ami2 > and Amu > 0), and tends to be 
brighter than the fourth arrival image (Amu < 0). Thus we 
expect that the first image to arrive is most likely the third 
brightest image; this provides support for our choice of the 
third brightest images for computing the magnification bias, 
because it assures the detection of the first lensed SN image 
to arrive. The appearance of this first SN image close to an 
early-type galaxy can act as an initial trigger, but two SNe 
are likely to be needed to confirm the system as a lens (and 
not a pair of SNe in the lens galaxy itself). The time delay 
between the first image and the subsequent images will still 
need to be measured, since the delay between images 2 and 
3 will often be very short. Our result indicates that typi- 
cal time delays for quad lenses range from ~ 1 day for the 
enshortest delay to ~ 1 month for the longest delay. 
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Figure 8. Predicted distributions of time delays and magnitude differences for various image pairs. Each panel shows values for 400 
image pairs randomly selected from our mock lens catalogue. Here we name images in order of arrival, i.e., image 1 arrives first, image 
2 next, etc. The upper panels show the distributions for lensed quasars, whereas the lower panels are for lensed SNe. 



5 COSMOLOGICAL PARAMETERS FROM 
TIME DELAYS IN LSST 

In this section, we use our mock LSST lens catalogue to 
explore one particular application of a future time-delay lens 
sample: the measurement of the Hubble constant, and also 
other cosmological parameters. 



5.1 Basic assumptions 

We follow the methodology proposed by lOguri ( 2007a]) to 
combine a number of time delay measurements. Th e tech- 
nique was adopted also bv lCoe fc Moustakasl (|2009T ) to dis- 
cuss future constraints from time delays. In brief, we intro- 
duce a reduced time delay defined by 



Atij 



2c A B 
1 + zi DiD s 



(19) 



where rj denotes the distance of image i from the center 
of the lens galaxy, and study its behavior as a function of 
the asymmetry Rij — \rj — ri\/\rj + rj| and the opening 



angle 8 i: 



[(xi ■ Xi)/(rirj)]. While the reduced time 



delay 3 becomes unity when the lens potential is exactly 



isothermal, external perturbations, non-isothermality, and 
substructures induce scatter in 3. 

We construct a subsample of lensed quasars and SNe for 
our statistical analysis from the mock lens catalogue. First, 
we restrict the image separation range to 0'.'8 < 8 < 2 / /5. The 
upper cut is meant to remove lens systems which are signifi- 
cantly affected by groups/clusters. We do not use lenses with 
< 0'.'8 because the detection and characterisation of lens- 
ing galaxies might be difficult for such small-separation lens 
systems. We also limit the lens redshift to z\ < 1.2 (see also 
Figure |5J, beyond which the 4000 A break of lensing galaxies 
comes in the z-band and therefore it becomes much more 
difficult to study the lensing galaxies without deep near- 
infrared images. Next, for simplicity we consider only the 
double image lenses that dominate the LSST lens sample. 
One of the reasons for this is our choice of the third bright- 
est images for computing magnification bias, which suggests 
that we may not be able to measure some of quad time 
delays without deeper follow-up monitoring. In either case, 
our forecast constraints are conservative in the sense that 
addition al quad lenses should on ly improve the constraints 
(see also lCoe fc Mou stakas 2009). Finally, we use only lenses 
with the asymmetry 0.15 < Rij < 0.8, in which the effects 
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of vari ous perturbations on S are modest and stable l|Oguril 
l2007af ). These cuts leave a sample of 1542 lenses (1476 lensed 
QSOs and 66 lensed SNe) for our forecast study. The sam- 
ple has average lens redshift z\ = 0.65 ± 0.27 (where we 
also give the standard deviation of the sample) and average 
source redshift z s = 2.29±0.86 (see also Figure[6]), which are 
reasonab ly similar to t he sim ple Gaussian distribut i ons as - 
sumed in lDobke etTafl (|2009l ') and lCoe fc Moustakasl (|2009h , 
although it should be kept in mind that these two redshift 
distributions are not independent but are naturally corre- 
lated. 

The error on S comes both from measu rement uncer- 
tainties and from mod el uncertainties (see lOguril |2007al : 
ICoe fc Moustakasll2009h . We assume measurement errors of 
cr(At) = 2.0 days and o"(r) = O'.'Ol. The assumed time de- 
lay measurement e rror i s the same as the one adopted by 
ICoe fc Moustakasl (|2009h . which was based on simulations 
by lEigenbrod et al.l (|2005l ) . We note that time delay mea- 



surements in the optical waveband are complicated by mi- 
crolensing for both lensed quasars (aslEigenbrod et al .11200 
showed) and lensed supernovae. iDobler fc Keetonl (|200' 
discuss the extent to which microlensing can bias time delay 
measurements in lensed supernovae. While we should be op- 
timistic about achieving 2-day precision given the wealth of 
additional imaging and spectroscopic data we can anticipate 
gathering on these precious systems, this is a topic requiring 
further research. 

Based on the Monte-Carlo simulations of l|Qguri2007al ). 
we adopt the model uncertainty of a(logH) = 0.08 which 
comes mostly from external shear and the scatter in the 
radial density slope. In this paper, we assume no error on 
the lens and source redshifts, which corresponds to the sit- 
uati on where all redshifts a re measured spectroscopically. 
See ICoe fc Mou stakas (2009) for the effect of photometric 
redshift errors on cosmological parameter estimates, which 
appears to be modest compared with lens potential uncer- 
tainties. 

As discussed in lOguril l|2007al ). the main source of sys- 
tematic error is the unknown mean logarithmic slope of 
the radial density profile of the lensing galaxies. The pa- 
rameter a — dlog<j)/dlogr for each lens is the effective 
slope of the lens potential, that includes the effect of the 
external convergence (a — 1 corresponds to an effectively 
isothermal potential, 0(x) oc r a ). We include the mean 
slope a as an additional nuisance parameter. We include 
the effect of changing a by adopting the scaling relation, 
Boc Afa oc (2-a) (e .g. jRefsdal fc Surdeil[l99l : IWitt et all 
l200d ; IWucknit3l2002l ). and assign it a Gaussian prior. We 
adopt the fiducial value of a pl ior(a) = 0.005, and marginalise 
over a to obtain constraints on cosmological parameters. 
We explicitly explore the impact of this assumption in Sec- 
tion [573] 



5.2 Fisher matrix analysis 

We obtain expected constraints on cosmological parameters 
using the Fisher matrix approach. We compute chi-squared 
as 



= E 



[log E(z h z s ) - log 5 mock ] z 
a(logH)2 



(20) 



Here cr(logH) is the quadratic sum of measurement and 
model uncertainties. The summation is over the 1542 mock 
lenses. Then the Fisher matrix for the parameter set {pi, 
P2, ■ ■ ■} is calculated as 

Here we adopt the standard parametrisation of the dark 
energy equation of state of the form: 

w(a) = w + (1 — a)iu a . (22) 

To complement the strong lensing data, we include the 
expected constraints from Cosmic Microwave Background 
(CMB) measurements and distance-redshift relation mea- 
surements from (unlensed) type-la supernovae. We compute 
the CMB Fisher matrix F™ B from upco ming Planck data 
using the simplified approach described in lMukheriee et"ail 
(2008). In the LSST era, we can expect this to be the 
basic prior PDF for dark energy studies l|Albrecht et al.l 
2006). For the SNe, we estimate expected constraints from a 
JDEM/SNAP survey assuming the redshift distribution o f 
SNe in th e final subsamp le given in lAldering et al.l |2004). 
Following iHu et all |2006), we add 300 local supernovae at 
2 = 0.05, and include both statistical and systematic errors 
in computing the Fisher matrix . In doing so, we include 
a nuisance parameter, M, which corresponds to the abso- 
lute magnitude of type-la supernovae, and marginalise over 
it. 

To summarise, assuming a flat universe, we consider the 
following 8 parameters: {wo, w a , S7 m , h, a, Qth 2 , n s , M}. 
Note that the time delays depend only on the first 5 of these 
parameters. Once the Fisher matrices are calculated, we ob- 
tain the parameter covariance matrix as C M „ = (F -1 ) M „, 
and display the corresponding marginalised constraints as 
ellipses containing 68% of the probability implied. 

5.3 Results 

We show the marginalised constraints on the cosmological 
parameters of interest in Figure [5] We find that the LSST 
strong lens time delays, combined with the Planck con- 
straints, provide marginalised la uncertainties (68% confi- 
dence intervals) of a(wo) ~ 0.15, cr(w a ) = 0.41, and cr(h) = 
0.017. If we include the JDEM/SNAP SNe constraints, we 
obtain a(w ) = 0.04, cr(w a ) = 0.17, and a(h) = 0.004. Al- 
though the constraints from time delays are less tight than 
those from the SNe, time delays still improve the constraint 
on the dark energy parameters because of the different de- 
generacy directions. This point is particularly clear in the 
wo-h plane, where the two constraints are almost perpen- 
dicular with each other. One way to view this is that future 
supernova surveys can be enhanced by strong lens time delay 
distance measurements. 

We also consider the case where we have a much weaker 
prior on the mean radial density slope a: if the prior PDF 
for this parameter is ten times broader than in our fiducial 
case, o"prior(a) = 0.05, the resulting constraints from time 
delays become much weaker (see Figure [9]) and hardly im- 
prove future constraints on dark energy from type-la SNe. 
This highlights the importance of accurate prior knowledge 
on lens potentials when deriving cosmological constraints 
from time delays. The slight rotation of the ellipse when the 
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Figure 9. Forecast marginalised constraints in the wo-Wn plane 
(upper panel) and the wo-h plane (lower panel). Ellipses show 
68% confidence regions for various combinations of expected con- 
straints from the Planck satellite (CMB), type-la supernovae in 
SNAP (SN), and time delays in LSST (At). Thin ellipses are con- 
straints from CMB+At, but assumes a much weaker mean slope 
prior of <T pr ; or (a) = 0.05 instead of the fiducial value of 0.005. 



prior on a is broadened is a result of the time delay dataset 
contributing less strongly to the joint fit - lack of knowl- 
edge of the mean density profile slope leads to lensing being 
downweighted as a cosmological probe. 

Indeed, the results presented here depend crucially on 
the width of the prior on the mean slope, o" pr ior(cv). In Fig- 
ure [TU] we show how the constraints on dark energy are af- 
fected by changing a pl i OI (a). We adopt the standard method 
of quantifying the constraints on the dark energy param- 
eters, employing the Dark Energy Task Force "figure of 
merit" (DE FOM), which is proport ional to the inverse o f 
the ellipse area in the uio-u> a plane |Albrecht et al.l l2006h . 
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Figur e 10. Dark energy figure of merit (FOM; e.g., 
lAlbrecht et alj|200 6f). defined as 1/ [cr(u) il )a(wp)] with w p being 
equation of state at the pivot rcdshift, as a function of the in- 
put (prior) uncertainty on the mean radial lens density profile 
slope a. The vertical dotted lines indicate our fiducial prior value, 



We find that time delays would not improve upon the con- 
straints from a JDEM/SNAP SNe survey if <7 prior (a) > 
0.01, which suggests that we need to determine the mean 
radial slope of galaxy density profile to better than ~ 0.01 
in order to use time delays in LSST as a future cosmologi- 
cal probe. In the limit of perfect prior information on a, we 
see that the LSST time delay sample could improve the DE 
FOM by a factor of 1.74, from 215 to 375. 



6 DISCUSSION 
6.1 Static surveys 

In this paper, we focused on cadenced surveys, mainly be- 
cause of the interest in strong lens time delays for use 
in measuring cosmological parameters, and the promising 
technique to find suc h lenses from their time variability 
IjKochanek et al.ll2006l ). As discussed in H3.2I however, some 
future wide-field optical imaging surveys will lack sufficient 
time domain sampling frequency to measure time delays. 
These surveys will still be useful for locating many lensed 
quasars, and possibly lensed SNe as well, but they will re- 
quire extensive follow-up monitoring observations to mea- 
sure the time delays. 

We note that some planned optical imaging surveys will 
be completely static, i.e., observing each field only once. Pre- 
vious lens surveys, however, have convincingly shown that 
we can still identify lensed quasar s by checking the colour 
and morphology of quasars (e.g., Ilnada et al.ll2008l ). High 
resolution, space-based surveys will be particularly efiec- 
tive here. For instan ce, the planned space mission Euclid 
jRefregier et al.ll2010h is planning to obtain high-resolution 
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imaging data for the entire southern extragalactic sky, which 
will be enormously useful for locating many lensed quasars 
from morphology. Once these quasar lens candidates have 
been identified, follow-up images can be taken to assess the 
variability of the candidates and to improve the efficiency 
of lens finding. Again, we will need additional monitoring 
data in order to measure time delays. Monitoring follow-up 
of lenses found in static surveys will, however, represent a 
significant observational undertaking. 



6.2 Follow-up observations 

Even for well-cadenced surveys such as LSST, we can con- 
duct follow-up monitoring to obtain more accurate time de- 
lays. In some cases the survey duration will not be enough to 
measure long time delays, and for these systems additional 
monitoring after the survey will be necessary. Particular at- 
tention should be paid to the follow-up of strongly lensed 
SNe because of their transient nature. Clearly it is impor- 
tant to identify lensed SNe as early as possible. In principle 
it will be possible to trigger follow-up monitoring by the 
appearance of the first lensed SN image, by checking the 
difference of photometric redshifts between the SN and its 
nearby galaxies, or from prior knowledge of the lensed host 
galaxies. Note that good image position and colour infor- 
mation for SN classification and redshift estimation will be 
required for this to work. Where these are not available, con- 
firmation of a lensed SN will likely require the appearance 
of the second image in a lensing configuration. 

In our Fisher matrix analysis, we assumed that all 
the redshifts are known. Spectroscopic redshifts for lensed 
quasars are relatively easy to obtain, thanks to prominent 
quasar emission lines. However, we can expect spectroscopy 
of the lensed SNe and the lensing galaxies to be generally 
much more difficult, because one has to detect absorption 
lines to measure the redshifts. Instead, we can adopt photo- 
metric redshifts, provided we can cleanly separate the lensed 
images from the lens galaxy light; if there is no bias inher- 
ited in the photometric redshift measurement, its effect is 
just to degrade the cosmological constraint s from the time 
delays slightly (see ICoe fc Mousta kas 2009). Spectroscopic 
follow-up observations are also helpful for the confirmation 
of strong lens systems, although the detection of time delays 
between the images does itself serve as a strong test of the 
lensing hypothesis. 



6.3 Prospect for constraining lens potentials 

One of our main findings f £|3.2[) is that accurate prior knowl- 
edge on lens potentials is crucial for cosmological constraints 
from time delays. Among others, the most important param- 
eter is the mean slope of the effective radial density profile of 
the lensing galaxy, a, including the effect of external conver- 
gence. Our requirement of a < 0.01 for the LSST time delay 
sample to be a useful addition to the cosmographic toolbox 
is in fact not far from the current s tatistical error of ~ 0.0 2 
from the Sloan Lens ACS Survey (jKoopmans et al.ll20Qgh . 
Care will need to be taken to ensure that the LSST lens 
galaxy sample is similarly well characterised. In addition, we 
expect that the external convergence may be well modelled 
by combining photometric information for the surrounding 



field with ray-tracing in iV-body simulations |Oguri et al. I 
120051 ; ISuvu et al.ll201(ih . 

With such a large number of strong lens systems, it 
may be possible to determine cosmological parameters and 
the mean radial density s lope simultaneously. For instance, 
iRusin fc Kochanekl (2005) constrained the mean radial den- 
sity profile by combining the Einstein radius measurements 
for lenses with different source redshifts. Moreover, we may 
also be able to use flux ratios between images to c o nstrain 
the r adial profile (e.g., iMortsell fc Sunessonl I2OO6I ; iMutkal 
2010), although we will need to take into account the ef- 
fects of dark halo substructures, differential dust extinction, 
and microlensing very carefully. We leave the exploration of 
these possibilities to future work. 



6.4 Comparison with individual lens modelling 

While we considered a statistical analysis for combining 
many time delay measurements, tight cosmological con- 
straints may also be obtained from the detailed modelling 
of individual lens systems (the so-called "golden lens" ap- 
proach). Indeed, th e feasibility of th is has been nicely 
demonstrated, by e .g TSuvu et all |2010l ) for B1608+656 and 
iFadelv et all (|2010r ) for Q0957+561. One advantage of our 
statistical approach is that it can handle various perturba- 
tions on lens potentials, such as substructures and higher- 
order external perturbations, more straightforwardly, since 
their effects average out when combining many lenses. Al- 
though the average does not always converge to zero (an 
example is external convergence), this poses no problems 
as long as the residual can be predicted reasonably accu- 
rately by theory. Indeed, the statistical framework provides a 
means, in principle, to infer the statistical properties of lens 
substructure, microlensing and dust simultaneously with the 
cosmology and lens profile parameters. 

For individual modelling we do not need to know the 
mean radial density profile very accurately, but instead need 
enough observational constraints to determine the radial 
profiles of each individual lensing galaxy. For this reason, we 
expect the quad lenses to be the focus of the individual mod- 
elling efforts, complenting the preference for doubles in the 
statistical approach. To summarise, we expect that individ- 
ual mass modelling will provide important complementary 
constraints on cosmological parameters, and help to assess 
the systematic errors that each methods involves. 



6.5 Group- or cluster-scale lenses 

Future surveys will discover not only galaxy-scale lenses but 
also wide-separation lenses produced by massive groups or 
clusters of galaxie s , such as those discovered in the SDSS 
(|lnada et al. 120031 . [2006). The time variability information 
will also be helpful in identifying these lenses as well, but 
the time delays between images will tend to be long, ~ 1 — 
10 years for massive clusters, suggesting that well-designed 
follow-up monitoring program may be necessary to measure 
their time delays. 
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7 CONCLUSIONS 

We have presented detailed calculations of the likely yields 
of several planned cadenced surveys, adopting realistic dis- 
tributions for the lens and source properties, and taking ac- 
count of the selection functions, including the magnifica- 
tion bias. We find that, for example, the LSST will discover 
~ 8000 lensed quasars (~ 3000 of which will have time delay 
measurements) and ~ 100 lensed SNe. Approximately one 
third of lensed SNe will be type-la. The lenses are dominated 
by double (two-image) lenses, with expected quad fractions 
of ~ 14% for lensed quasars and ~ 30% for lensed SNe. 
We have also produced mock catalogues of lenses, which are 
useful for probing strong lensing selection effects and the 
feasibility of various science projects. 

We have used a mock catalogue of ~ 1500 well-observed 
double lenses in LSST to derive expected cosmological con- 
straints. Specifically, we derived precisions on the Hubble 
constant and the dark energy equation of state parameters 
from this sample of time delay measurements, assuming pri- 
ors from Planck. The resulting predicted marginalised 68% 
confidence intervals are a(wo) — 0.15, cr(w a ) — 0.41, and 
u(h) = 0.017, implying that time delays can improve con- 
straints from a JDEM supernova type-la sample. However, 
this result holds only if we have accurate prior knowledge of 
the lens population's mean effective density profile: we find 
that the prior on the mean radial density slope has to be 
q < 0.01 in order for time delays from the LSST strong lens 
to be a useful future cosmological probe. 
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